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AERONAUTIC SYMBOLS 
1. FUNDAMENTAL AND DERIVED UNITS 
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2. 
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Length 
Time_- 
Force_. 

Power _ 
Speed . 



Symbol 



I 
t 

F 



P 
V 



Metric 



Unit 



meter 
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weight of 1 kilogram _ 



horsepower (metric), 
("kilometers per hour. 
\meters per second- _ 



Abbrevia- 
tion 



m 
s 



k.p.h. 
m.p.s. 



English 



Unit 



foot (or mile) 

second (or hour)., 
weight of 1 pound 

horsepower 

miles per hour 

feet per second 



Abbrevia- 
tion 



ft. (or mi.) 
sec. (or hr.) 
lb. 



hp. 

m.p.h. 
f.p.s. 



Weight = mg 

Standard acceleration of 
m/s 2 or 32.1740 ft./sec. 2 

Mass = — 
9 

Moment of inertia = mk 2 . 

radius of gyration k by proper subscript.) 
Coefficient of viscosity 



2. GENERAL SYMBOLS 

v, Kinematic viscosity 



gravity = 9.80665 



(Indicate axis of 



p, Density (mass per unit volume) 

Standard density of dry air, 0.12497 kg-m~ 4 -s 2 at 

15° C. and 760 mm; or 0.002378 lb.-ft." 4 sec. 2 
Specific weight of "standard" air, 1.2255 kg/m 3 or 

0.07651 lb./cu.ft. 



Area 

Area of wing 
Gap 
Span 
Chord 

Aspect ratio 

True air speed 

Dynamic pressure : 



3. AERODYNAMIC SYMBOLS 

i m Angle of setting of wings (relative to thrust 



to 

n, 

VI 



\pV 2 



Lift, absolute coefficient C L = ^ 



D 0 



Drag, absolute coefficient ^ — ^ 

Profile drag, absolute coefficient C Do = ^ 

Induced drag, absolute coefficient G Di ^~^ 

Parasite drag, absolute coefficient = 

Cross-wind force, absolute coefficient (7 c a 
Resultant force 



C_ 



C p , 

*, 

Ot 0 , 
<Xa, 
7i 



line) 

Angle of stabilizer setting (relative to thrust 
line) 

Resultant moment 
Resultant angular velocity 

Reynolds Number, where I is a linear dimension 
(e.g., for a model airfoil 3 in. chord, 100 
m.p.h. normal pressure at 15° C, the cor- 
responding number is 234,000; or for a model 
of 10 cm chord, 40 m.p.s. the corresponding 
number is 274,000) 

Center-of-pressure coefficient (ratio of distance 
of c.p. from leading edge to chord length) 

Angle of attack 

Angle of downwash 

Angle of attack, infinite aspect ratio 

Angle of attack, induced 

Angle of attack, absolute (measured from zero- 
lift position) 
Flight-path angle 
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SUMMARY 

The speed distribution in a laminar boundary layer on 
the surface of an elliptic cylinder y of major and minor 
axes 11.78 and 3.98 inches, respectively, has been deter- 
mined by means of a hot-wire anemometer. The direction 
of the impinging air stream was parallel to the major axis. 
Sjwcial attention was given to the speed distribution in 
the region of separation and to the exact location of the 
point of separation. An approximate method, developed 
by K. Pohlhausen for computing the speed distribution, 
the thickness of the layer, and the point of separation, is 
described in detail; and speed-distribution curves calcu- 
lated by this method are presented for comparison with 
experiment. Good agreement is obtained along the for- 
ward part of the cylinder, but Pohlhausen' s method fails 
shortly before the separation point is reached and conse- 
quently cannot be used to locate this point. 

The work was carried out at the National Bureau of 
Standards with the cooperation and financial assistance 
of the National Advisory Committee for Aeronautics. 

INTRODUCTION 

PrandtFs postulate, that the viscosity of the fluid 
plays an important role only in a thin " boundary 
layer" on the surface of a body in a moving fluid, has 
been the key to the understanding of many puzzling 
phenomena occurring in the flow of viscous fluids. 
Subsequent works, both experimental and theoretical, 
have verified the existence of a " boundary layer" and 
have given conclusive evidence that the radical differ- 
ence between the behavior of an ideal fluid and that of 
a viscous fluid is due to the presence of a boundary 
layer and the wake resulting from its separation from 
the surface of the body. 

In the region of potential flow outside the boundary 
layer and wake the effect of viscosity is negligible; 
hence the mathematical theory pertaining to the flow 
of a nonviscous fluid may be applied, a procedure 
which obviously finds greatest application in the treat- 
ment of flow about body forms on which separation is 
delayed and for which the wake is small. While the 
flow within the boundary layer must be treated as vis- 
cous, much of the difficulty was avoided by Prandtl, 
who derived a general equation of motion on the | 



assumption that the boundary layer is thin compared 
to the dimensions of the body. Thus, excluding the 
wake, it might appear possible to work with two re- 
gions, one being the boundary layer and the other the 
region outside, each governed by its own set of laws. 
A more careful consideration shows, however, that 
each has an effect upon the other, and that no satisfac- 
tory solution of the problem as a whole can be obtained 
without treating the two simultaneously. At present 
this is impossible. 

It has been found that the velocity distribution in 
the boundary layer, its thickness, and its tendency to 
separate from the surface of the body are governed 
almost entirely by the velocity distribution in the 
region of potential flow outside the layer. The pro- 
cedure adopted, therefore, is to measure the velocity 
distribution outside the layer and to use this as a basis 
for calculating the flow within the layer itself. This 
procedure takes the place of the familiar one used in 
potential theory, that of using the shape of the body as 
the basis. A treatment of the former type, employing 
the Karman integral equation, has been developed 
by K. Pohlhausen (reference 1). Unfortunately this 
treatment is approximate and experimental data are 
necessary to test the validity of the assumptions made. 

The flow in the boundary layer described by 
Prandtl's equation is termed "laminar flow." It is 
well known that as the layer thickens this type of flow 
breaks down and is replaced by turbulent flow. If it 
be desired to maintain laminar flow, transition to the 
turbulent state may always be avoided by a sufficient 
lowering of the Eeynolds Number. 

The present paper describes a study of the flow in a 
laminar boundary layer formed in the two-dimensional 
flow around an elliptic cylinder oriented with its major 
axis parallel to the air stream of a wind tunnel. Meas- 
urements of the speed distribution in the layer have been 
made by means of a hot-wire anemometer, and extend 
from the origin of the layer to and beyond the region 
where it separates from the surface of the cylinder. 
The Reynolds Number has been so chosen that the flow 
in the boundary layer does not become turbulent before 
separation. The purpose of the study was to test 
Pohlhausen's approximate solution of the Karman 
integral equation, especially in the region of separation. 

l 
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While a treatment of boundary-layer theory may 
readily be found in the literature (for example, in 
reference 1, which includes also Pohlhausen's solution) 
a brief discussion is included for conveniently judging 
the validity of the approximations as applied to the 
problem at hand. 

L DEFINITIONS OF SYMBOLS 

We shall adopt a curvilinear system of coordinates in 
which x denotes distance along the curved surface of the 
elliptic cylinder measured from the forward stagnation 
point and y distance normal to and measured from the 
surface. The radius of curvature of the surface, and 
hence of the x axis, will be denoted by r. Since y will 
always be small, we shall neglect its curvature. The x 
and y components of velocity in the boundary layer 
will be denoted by u and v, respectively, and the x com- 
ponent of velocity just outside the boundary layer, by 
U. The pressure of the air at any point in the layer 
will be denoted by p. 

We shall find it convenient at the outset to express 
distances, speeds, and pressure in terms of correspond- 
ing reference quantities as the units, denoted respec- 
tively, by X, U 0 and P 0 . The minor axis of the ellipse 
(L=3.98 inches) is selected as the reference length. 
The speed of the undisturbed stream of the wind tunnel 
is chosen as the reference speed; that is, the air speed 
that would prevail in the empty tunnel at the position 
to be occupied by the leading edge (forward stagnation 
point) of the elliptic cylinder (U 0 , about 11.5 feet per 
second). The reference pressure (really a pressure dif- 
ference) conveniently follows as the pressure rise when 
air of initial speed U Q is brought to rest by impact, 

I e ., P 0 =% p U 0 2 , where p is the air density. The 

Reynolds Number R is defined as where v is the 

kinematic viscosity of the air. 

The foregoing quantities are made dimensionless as 
follows: 

x y r_ u_ v_ JL- V 

T L' If Uo Uo P 0 ~HpW 
In order to avoid writing quotients, the symbols x, y, 
r, u, v, U, and p will be used to denote dimensionless 
quantities in the remainder of the paper, the division 
by the appropriate reference quantity being inferred. 
The symbols with their new meaning are summarized 
as follows: 

L, Length of minor axis of elliptic cylinder. 

x, Distance along the curved surface of the elliptic 
cylinder at right angles to the axis of the cyl- 
inder and measured from the forward stagna- 
tion point, divided by L. 

y } Distance normal to and measured from the 
surface divided by L. 

r, Kadius of curvature of the surface divided by 
L. 



u, x component of velocity in the boundary layer 
divided by U 0 . 

U ) x component of velocity in the region of poten- 
tial flow just outside the boundary layer 
divided by U 0 . 

v, y component of velocity in the boundary layer 
divided by U 0 . 

p, Pressure of the air divided by P 0 . 
For describing the separated boundaiy layer x b and 

?/ b are introduced. 

x b) Distance divided by i, the same as x up to 
the point of separation, but parallel to the 
major axis of the ellipse beyond this point. 
y b) Distance divided by L perpendicular to x b , 
with origin on the major axis of the ellipse or 
the major axis produced. 
(x b) y b ) is any point aft of the separation point. 
Figure 13 may help the reader to visualize some of 
the foregoing quantities. 

II. BOUNDARY LAYER THEORY 

1. GENERAL EQUATION 

The experimental conditions, a steady two-dimension- 
al flow of air at speeds so low that compressibility may 
be neglected, makes possible an immediate simplifica- 
tion of the Navier-Stokes equations (reference 2). 
Their simplified and appropriate form, given by Fuchs- 
Hopf (reference 3) in our particular curvilinear co- 
ordinate system and written here nondimensionally, is 



du du uv l/d 2 u . d 2 ' 
u dx + % r ~R\dx 2± dy 2 

_1 
2 



d?u_l 3^_2 dv_u , v_ dr\ 
dy 2 ~ r dy r dx r 2 r 2 dxj 



dp 
dx 



(1) 



dv . dv , u 2 _ l/d 2 v d 2 v 2 du I dv_u dr_v\ 
u ^ + % + 7~R\dx 2+ dy 2± r dx^ r dy r 2 dx r 2 ) 



(2) 



(3) 



_1 dp 
2 dy 

These together with the equation of continuity 

du . dv^v^Q 
dx dy r 

are the equations of motion of a viscous fluid from 
which a mathematical treatment of boundary-layer 
flow may begin. 

A boundary-layer equation was derived by Prandtl 
(reference 4) from equations in rectangular coordinates 
corresponding to (1), (2), and (3) by neglecting terms 
that become unimportant when the treatment is ap- 
plied to a thin layer of thickness 5 (5 is also expressed 
with L as the unit). The terms that may be neglected 
were found by a consideration of relative orders of 
magnitude. Although the quantities retained arc on 
the average, of a higher order of magnitude than those 
neglected, in certain regions some of the former assume 
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very small values and even change sign. In these 
regions some of the neglected quantities may be of a 
higher order than some of those retained but, from the 
form of the equations, it does not appear that this will 
introduce any significant additional uncertainty. 
The ultimate test of the accuracy of the approxima- 
tion formulas is the agreement between calculation and 
observation. 

In the following discussion we shall use the expres- 
sion "of the order of" (designated by the symbol 0) to 
mean the maximum order attained by a given quantity. 
For example, when u is said to be 0 (1), u may reach 
values ±1 or even ±10, but it may at times have 
much smaller numerical values, even zero. A factor 
of 10, or possibly greater, is not to be considered as 
sufficient to change the order of magnitude but a 
factor of, let us say, 100 places the quantity in a higher 
older. 

Experimentally we find sufficient justification for 
writing 

. TT du d 2 u , dU n /in 

x ' r ' u > u > Fx a? and & are0 (1) 



du . 



dy 



is 0 



(i) 

W 15 °W 

From Bernoulli's theorem (p+ U 2 = constant, bear- 
ing in mind that this equation is in dimensionless 
form) 

t " o CO 

dv v 

Equation (3) indicates that - is 0 (1). Since 

dv 

v can safely be assumed to be 0(5), must be 0(1). 

Consideration of the order of mangitude of v indicates 
that it is a reasonable assumption to assign values 

00) to g£ and and a value 0^- J to ^ 

Equations (1) and (2) are now rewritten with orders 
of magnitude replacing the terms themselves, except 

for R and which are as yet undetermined. 

0(l)+0(l)-0(5)=g[0(l) + o(i)-o(i)-0(5)-0(l) 

+ 0(5)]-0(l) ' (la) 



O(«) + O(«) + O(l)=g[O(5) + O0)+O(l) + O(l)-O(l) 



(2a) 



i r for the elliptic cylinder goes from 0.17 at the leading edge to 4.4 at the thickest 
section. 



In the boundary layer the inertia! and I fictional forces 
are assumed to be of nearly equal importance; hence 
the inertia! and frictional terms in equations (1) and 
(2) must be of the same order of magnitude. In order 
to satisfy this requirement, 5 must be such that R in 

equation (la) is O^Jjy; and, to satisfy equation (2a) 

at the same time, must be 0 (1). 

The assumption of a thin layer means thai 8 is small 
compared to 1, and hence that terms 0 (5) or smaller 
may be neglected. When this is done, equation (1) 
becomes 



du 



dx~*~ dy~~R dy 2 2 dx 



(4) 



and equation (2), while of no importance in the dis- 
cussion to foll( 
0 (1), becomes 



cussion to follow other than to indicate that is 

dy 



2 dy~ r 



Equation (3) reduces to 



du^dv q 

dx dy~ 



(5) 



Equation (4) is the boundary-layer equation as orig- 
inally derived by Prandtl. 

Two additional results have come out of the preced- 
ing discussion. The first, that is 0 (R) } not only 

indicates one experimental requirement for the validity 
of equation (4) but shows us, to order of magnitude at 

least, that 5 varies as when equation (4) is valid. 

By momentum considerations the law is actually found 

to be 6 = K^j£> where K is an unknown constant of 

0 (1). The second, that g| is 0 (1), shows that the 

change in pressure across the boundary layer from 
?/=0 to y= 5 is 0 (5). The significance of the latter 
deduction is that a pressure measured at some point- 
on the surface of the cylinder may be regarded as the 
pressure in the boundary layer and in the potential 
flow just beyond it, anywhere along the normal to the 
surface from that point, and hence may be used to 
compute the speed in the potential flow by Bernoulli's 
equation. 

At the surface u and v are zero. As y increases, 
U asymptotically approaches £7, the speed in potential 
flow outside the boundary layer. U is usually a 
function of x; and just beyond the boundary layer 
may be regarded as a function of x alone, its change 
with ?/ in the potential flow being of a lower order of 
magnitude than within the boundary layer. 
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In potential flow we have by Bernoulli's theorem 

P+U 2 =p 0 +U 0 2 

where p and p 0 are the dimensionless pressure for the 
points with the dimensionless velocities U and Uq, 
respectively. We shall choose for U 0 the speed of the 
undisturbed stream; and when we do so, we note, upon 
reference to section I, that its value is unity. We 
observe also that p 0 is a constant. Replacing U 0 by 
unity and solving for i7, we obtain 

U 2 =l-(p-po) (6) 

Equation (6) is a relation between the dimensionless 
pressure and the dimensionless speed at the outer 
boundary of the layer. Differentiating equation (6) 

with respect to x and denoting ^ by U' yields 



J*=2VTJ' 
ax 



(7) 



Neglecting a pressure difference 0 (5), the pressure at 
the surface and its variation with x are related to U 
and U' by equations (6) and (7). 

Using equations (4), (5), and (7) an integral equation, 
namely, the "Karman integral equation", may he 
developed (reference 1). It is expressed as follows: 



>.U'J\u-u)dy- Uj x j \u-u)du 



+ 



(8) 



y=0 



_ POULHAUSEN'S APPROXIMATE SOLUTION OF THE KARMAN 
INTEGRAL EQUATION 

Equation (8) can be solved only when some relation 
between v and y is assumed. Pohlhausen chose 



-Ay+By 2 +Cif+Dy A 



GO 



where A, />, < ', and I) are determined by (he boundary 
conditions at the surface and al the outer limit of the 



o.s[- 



layer. A parameter 8 is introduced such that at 
// = 5, u=U. Having defined 5 in this way, we may 

I) f 

regard it as the boundary-layer thickness. 



Since -5- 
oy 



du 



-9072 + l()70.4X 



is of a lower order than -^y we may write at y=8 } 
du d U 

=—=0. The supposition has been made that the 
dy dy 

thickness has a definite finite value; whereas in the 
physical case u approaches U asymptotically, there 
being no sharp division between the boundary layer 
and the region of potential flow. The conditions at 
5 may be made to simulate the asymptotic case more 
closely by imposing the condition that the curve of 
equation (9) shall reach the line u=U with zero curva- 
ture. As outer boundary conditions we have therefore 
TT du A d 2 u ft , . 
u=Uj &y ==0) W 2 afc ^ = 5 
As inner boundary conditions we have 

d 2 u R dp TJT1 , jy n 
dtf = 2dt = - UU R * ty = ° 

by equations (4) and (7), since u=0 and 6=Q. 

When these four conditions are applied, equation (9) 
becomes 

KK«H-^iXrK>-*X9> 

where X= U' 8 2 H. Equation (10) represents a family 
of curves, with 5 as parameter, satisfying the inner and 
outer boundary conditions. 

In keeping with the foregoing conditions, the upper 
limit of the integrals in equation (8) may be written 8 
instead of 00 since U—u is zero when y is greater than 8. 
When this is done and u as given by equation (10) is 
introduced into equation (8), the integrations with 
respect to ?/ may be performed. Introducing, tor 
convenience, \=*U* z and z=R 8' 2 and making the 
differentiations with respect to x indicated in (8), 
the variables being U, I 1 ', and 2, we obtain finally 



4 + 4.8- 



(11) 



dx~ [/(-213.12 + 5.76X+X 2 ) 

and E7" are known functions of x, equa- 



When U, I 

(ion (1 1 ) may be solved for z. The method of obtain- 
ing the solution is given in section III. For example, 
at some particular value of x, say x ti Zj is found by 
solving (11). From z=R 8 2 and \=U' z the partic- 



ular VJi 



lines 0 



and \i=U\ z u where I '', is the 



value of U' at X\ } may be found and substituted in 
(10). The result is 



u=U l 

where l7 2 is the value of U at ^. Equation (12) is a 
relationship between u and y^R at the point 

Since Xi, Z\ and 5, are constants and y^/R is V^f*' 
the (y-yjRy u ) curve is the same for all values of 



(12) 

the Reynolds Number as long as the assumptions 
underlying the boundary-layer equations are valid. 
Since the true curve can be represented only approxi- 
mately by a fourth-degree equation, equation (12) 
can lead to only an approximate solution. 
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Separation of the flow from the surface occurs when 
the fluid in the inner part of the boundary layer has 
lost sufficient kinetic energy through friction to fail 
to advance farther into a region of rising pressure. 
Back of the position where the fluid comes to rest, 
the higher pressure downstream causes reversed flow 
and the consequent accumulation of fluid crowds the 
main flow away from the surface. The point where 
the fluid comes to rest, just back of which the flow 
reverses in direction, is the separation point. Analy- 
tically, separation will occur when the initial slope of 



and is 54 inches in length. It was mounted verti- 
cally between the upper and lower faces of a 54-inch 
octagonal wind tunnel, thus extending completely 
across the tunnel from one face to the opposite face. 
The surface was made smooth by alternate varnishing 
and sandpapering and finally by polishing with crocus 
paper. Around one side 2 l /2 inches below the hori- 
zontal mid-section of the tunnel, 16 pressure orifices 
were inserted for obtaining the pressure distribution 
over the surface. The connections were copper tubes 
of 0.04-inch inside diameter extending from the outer 




Figure 1. — Photograph showing elliptic cylinder and hot-wire anemometer. The microscope as it was set for viewing the wire is also shown. 



the curve representing equation (10) changes from a 
positive to a negative value. This occurs when the 
initial slope is zero; that is, when the first term of 
equation (10) is zero. The criterion for separation 
resulting from Pohlhausen's solution is then 
\=U'8 2 R=-i2. 

III. EXPERIMENTAL STUDY OF BOUNDARY LAYER 

1. DESCRIPTION OF ELLIPTIC CYLINDER AND HOT-WIRE 
ANEMOMETER 

The elliptic cylinder is of wooden construction, with 
major and minor axes of 11.78 inches and 3.98 inches, 



surface to the hollow interior of the cylinder and thence 
to the outside of the tunnel where connection could 
be made to a manometer. The number and distance 
of each orifice measured along the surface from the 
stagnation point are given in the first and second 
columns of table I. 

The hot-wire anemometer may be described best by 
reference to the photograph (fig. 1). The essential 
part of the instrument is a platinum wire, 0.002 inch 
in diameter and 3 inches in length, on the ends of the 
flexible supporting prongs. The heat loss from this 
wire when heated by a constant electric current serves 
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as an indication of the air speed across it. By means 
of the three micrometer screws in the small lathe head 
on the bracket at the rear of the cylinder, the wire 
could be moved in any direction through a distance 
of about 2 inches. A setting to any position around 
the surface could be made by readjusting the upper 
horizontal part of the support. 

The microscope and lamp shown in figure 1 were 
used to determine the distance between the wire and 
the surface, for a given micrometer reading, by Mew - 
ing the wire and its reflection in the surface. They 
were removed when speed measure incuts were being 
made. 

2. PRESSURE DISTRIBUTION 

An inclined benzol manometer, one side of which was 
open to the still air of the tunnel room and the other 
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FIGURE 2.— Pressure distribution around one side of elliptic cylinder as determined 
at 16 pressure orifices. 

side connected to a calibrated pressure orifice in the 
tunnel wall 10 feet ahead of the cylinder, was used in 
measuring air speed in the tunnel. The movement of 
the liquid in the manometer was slight, since the 
speed maintained was low (about 11.5 feet per second i, 
and had to be measured by a traveling microscope 
provided with a micrometer screw. The difference 
in pressure between the wall of the tunnel and each 
of 16 positions about the cylinder was obtained by 
connecting separately the pressure orifices in the cyl- 
inder surface to the side of the manometer previously 
open. The prevailing drop in static pressure between 
the tunnel-wall orifice and the leading edge of the 
cylinder was subtracted from these differences to give 
p—jh, from which U was computed by means of 
equation (6). 

After a setting of the cylinder with its major axis 
approximately parallel to ihe undisturbed air stream, 



finer adjustment of the orientation was made until the 
angle of attack was zero as indicated by comparison 
between the observed pressure distribution over the 
forward part of the cylinder and that computed from 
nonviscous fluid theory for zero angle of attack (refer- 
ence 5). Exact agreement with theory could not be 
obtained because of the influence of the boundary layer, 
the blocking of the tunnel stream by the cylinder, 
and the pressure drop in the tunnel. For making 
boundary-layer calculations, however, only the actual 
pressure distribution need be known. It may be re- 
marked that the purpose was to study not the boundary 
layer on the surface of an elliptic cylinder, but rather 
the boundary layer on a surface with a known pressure 
distribution. 

The effect of the hot-wire anemometer and support 
upon the distribution was determined for various set- 
tings about the surface with the anemometer lowered 
until the prongs spanned the orifices. Since the effect 

7f 
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Figure i— Theoretical and observed speed distributions at stagnation point . The 
theoretical curve was obtained from ideal fluid theory for the elliptic cylinder at 
zero angle of attack. 

was small in all cases, especially ahead of the prongs, 
it was neglected. Pressure distributions were measured 
from time to time during the work as a safeguard against 
accidental shifting of the cylinder. Good agreement 
was always obtained. The final results for anemometer 
and support absent are given in column 3 of table I 
and in figure 2. In figure 2 the theoretical pressure 
distribution for flow of a nonviscous fluid is also given. 

3. SPEED DISTRIBUTION IN THE BOUNDARY LAYER 

The relationship between the rate of cooling of a 
heated wire and the speed of the air across it is 

H=A+BjS 

where II is the heat loss per degree difference in tem- 
perature between the wire and the air (in the present 
experiment, expressed in watts per degree Centigrade), 
Sis the cross-wire air speed, and A and B are functions 
of the particular wire and the air density. .1 and B 
are determined by measuring H at known air speeds. 
When // is plotted as ordinate and 7$ as abscissa, 
the resulting calibration curve is linear. Since A and 
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B were found to change gradually over a period of 
time, or suddenly if the wire received rough treatment, 
it was necessary to calibrate the wire before and after 
each day's work. Since the hot wire responds equally 
well to flow in any direction, only the resultant speed 
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Figure 4.— Speed distributions in the boundary layer at three positions as desig- 
nated by the values of x. Note the better approach to the origin when no correc- 
tion is made for heat loss to the surface. 

can be measured with it. For this reason u and U 
are termed "speeds" rather than "velocities" when 
speaking of experimentally determined values. Before 
the boundary layer separates from the surface the 
measured speeds are, with good approximation, the 
x components of velocity. Wherever the flow direction 
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is obviously not that of x, the superscript * will be 
added, for example, u*. 

Before measurements of speed distribution were 
made, the effect of the cylinder surface on the heat 
loss from the hot wire in still air was determined. 
At 2 millimeters from the surface the cooling effect 
became perceptible. Within this distance the heat 
loss to the surface (expressed as H—II™ , where Hcq is 



the heat loss in still air at a great distance from the 
surface and H is the heat loss near the surface) was 
determined as a function of the distance y and the 
temperature difference between the wire and the 
surroundings. It was recognized at the time that the 
heat loss to the surface is also a function of the air 
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Figure 6.— Speed distributions in the boundary layer. 

speed at the wire and the gradient in speed between 
the wire and the surface. Some results of Piercy 
and Richardson (reference 6) strengthen this view. 

means were available, however, for determining 
this relationship, nor were the results of Piercy and 
Richardson applicable. Just how and w r hen a cor- 
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Figure 7.— Speed distributions in the boundary layer. Note the better approach 
to the origin when a correction is made for heat loss to the surface. 

rection should be applied to eliminate heat loss to the 
surface is therefore uncertain. When the smooth 
approach of the speed to zero at the surface is used 
as a criterion, the present results indicate that a 
correction varying as some inverse function of the 
speed, as well as with distance and temperature 
difference, should be applied. 



10 



8 



REPORT NATIONAL ADVISORY COMMITTEE FOR AERONAUTICS 



Traverses through the boundary layer were made 
at suitably chosen positions on the surface. At each 
position measurements of the speed were made from 
well outside the boundary layer to within 0.2 milli- 
meter of the surface. The intervals moved were read 
from the scale attached to that one of the three 
micrometer screws which was employed at the time 
in moving the wire. The final distance y was deter- 
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Figure 8— Speed distributions in the boundary layer. Note the better approach 
to zero when a correction is made for heat loss to surface. 

mined by counting divisions on the eyepiece scale of 
the microscope between the center of the wire and 
its image in the surface. Each time the microscope 
was set into position for measuring this final distance, 
the eyepiece scale was calibrated by comparing the 
change in the number of scale divisions between the 
wire and its image with the distance through which 
the wire was moved by the screw. It is estimated 
that distances so determined have an accuracy of 
± 0.04 millimeter. 

The results are given in figures 3 to 12, inclusive. 
By means of the (y^lR, u) scheme of representation 
the curves are independent of Reynolds Number. 
(See sec. II, 2.) Reference to the x and x h scales of 
figure 13 will aid the reader to visualize the position 
of the various distributions. It will be observed that 
two sets of results are given, those to which a correc- 
tion for heat loss to the surface as determined in still 
air has been applied and those for which it has been 
omitted entirely. Allowing for an uncertainty in 
y^R of ±0.06 (±0.04 mm), the manner of approach 
to zero will aid in judging where the correction should 
be applied and where it should not. In figures 4 
and 5, because of the high speeds at which the lower 
points fall, the correction seems to be unnecessary; 
whereas in figures 6 to 12, with somewhat lower 
speeds at the lower points, a correction seems to be 
justified. 

In figure 3 a comparison is made between the 
measured speed distribution at the stagnation point 



and that computed on the basis of an ideal fluid 
(reference 5). Although exact agreement is not to be 
expected because of the tunnel blocking and the effect 
of the wake, such a comparison serves as an approxi- 
mate check on the hot-wire measurements. 

A clear and precise indication of the position of 
separation was obtained by allowing smoke to enter 
the stream slowly through a pressure orifice just back 
of the separation point. The smoke moved slowly 
forward (indicating reversed flow) to a definite stop- 
ping point. When another orifice farther back of 
separation was employed, the smoke from it moved 
forward to the same limiting position. The separa- 
tion point was thus determined at x= 1.99 ±0.02. 
The experimental curve of figure 8 shows the speed 
distribution just before separation and that of figure 9 
just after. That no marked change appears may be 
attributed to the fact that the hot-wire anemometer 
responds equally well to flow in any direction. This 
uncertainty in direction is indicated by the use of a 
special symbol w*, it being understood that the un- 
certainty exists only on the lower horizontal part of 
the curves representing flow after separation. Before 
separation, with u of order 1 and v of order 6 = 0.05, 
the assumption that the measured speed is parallel 
to x is quite good. After separation not only has the 
velocity reversed in direction near the surface but, in 
addition, there are probably large normal components 
where the reversed flow curves outward. If we were 
to speak of velocities in connection with the experi- 
mental curve of figure 9, we would say that the 



































































































s 




































P 

s 


r 




































































/ 


/ 










— Pohlhausen's 
















i 




/ 












approx/more 
Solution. 
Observed (of 


















/ 

/ 

/ 










-( 


)-£ 


fer 










/ 


/ 

/ 
















correction mode 








7 


i 


o' 
















for heot los 
to surface, 
^orrection.o 


5 












/ 

/ 

/ 














) 


< C 


5 


1 

0 








o 


















aererminea in 
still air. ooo/iec 


V 






-7/7, 


ner 


■ fin 


lit i 


Of 6 


•ep 


ore 


it- 


x^2.029. { R 


=2: 


16 C 


/ 1 




ea Dounaary /oyer 

\ 



















0 123-456789 

Figure 9.— Speed distributions in the boundary layer after separation. The * 
indicates an uncertainty in the direction of the speed in the wake. 



vector representing the velocity turns from a direc- 
tion parallel to the surface pointing upstream through 
a little less than 180° pointing downstream in approxi- 
mately the direction of x b , as we follow the curve 
from the origin outward. 

Figures 9, 10, 11, and 12 show traverses of the sepa- 
rated boundary layer and wake. It is apparent that no 
abrupt change in speed distribution in the layer accom- 



AIR FLOW IN A SEPARATING LAMINAR BOUNDARY LAYER 



9 



panies separation. It must be remembered, however, 
that u* is resultant speed and that approximately the 
same conditions as to direction obtain in figures 10, 11, 
and 12 as in figure 9. Visual observations with smoke 
showed that the air motion in the wake fluctuated 
greatly and that only near the separation point was 
the reversed flow at all regular. The hot wire shows 
this persistence of air motion in the wake. The results 
ol >( .-lined here are much like those obtained by Linke 
in his study of the laminar boundary layer separated 
from a circular cylinder (reference 7). In 
figure 13 the boundary layer is drawn to the 
same scale as the elliptic section. The thick- 
ness of the layer both on and separated from 
the surface is taken as the width of the region 
in which the speed is changing; that is, the 
length of abscissa between the points of the 
(yjR, u) curves where the tangents become 
horizontal. Owing to the asymptotic char- 
acter of the limits, 8 is somewhat arbitrary; 
but, since the same procedure was used 
throughout, figure 13 gives a picture of one 
part of the boundary layer as related to any 
other and shows the changes that occur. 

IV. CALCULATION OF SPEED DISTRIBUTION 
BY POHLHAUSEN'S APPROXIMATE METHOD 

The information required for an evaluation 
of U, U', and -jpr as functions of x may all be 



small correction. For values of x greater than 1.6, where 
the theoretical and experimental curves bear no resem- 
blance to one another, a graphical determination of 
the slope of the experimental (x } p—jh) curve was made 
directly; and for values of x between 1.1 and 1.6 both 
the direct and difference methods were used. A curve 

of (x, was then plotted and a similar procedure 

d 2 v 

followed in determining The final values are 
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Figure 11.— Speed distributions in the boundary layer after separation. 

uncertainty in the direction of the speed in the wake. 
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The * indicates an 



obtained from the curve of pressure distribution and 
its first and second derivatives. The derivatives were 
obtained by the following procedure. From x=0 to 




given in columns 4 and 5 of table I. By the formulas 

(a), (6), and (c) at the foot of table I, U, U' and 

were calculated. Formula (a) is the same as 
equation (6) and (b) and (c) follow from (a) 
by differentiation and combination . The val- 



ues of U, U' 9 and 



UU" 



are given in columns 



U' 2 

6, 7, and S of (able I. 

Equation (11) was solved by the isocline 
method. This method involves computing 

^ with assumed values of z for various val- 
ues of x and the accompanying values of U, 



V' and 



^ytz ' Points of known ^ are then 



5 6 

Figure 10.— Speed distributions in the boundary layer after separation. 

uncertainty in the direction of the speed in the wake. 



The * indicates an 



1.1, where the theoretical and experimental pressure 
distributions agree rather closely, a curve of differences 
was plotted with x as abscissa. The slopes of the dif- 
erence curve, obtained graphically, were then applied 
to the theoretical slopes as a correction. In this way 
errors in determining the true slope entered only into the 



located on an (x, z) diagram and a short line 
drawn through each to indicate the slope at 
the point. The solution curve is then drawn 
in, passing through the various points with 
the slope indicated at each. (See fig. 14.) 
A procedure more convenient in certain parts 
of the diagram is to draw lines of constant 
slope (isoclines) and then draw the solution curve 
crossing the isoclines with the slope marked on them. 

The particular solution curve required is the one 
that satisfies the boundary conditions at x=0, the 

dz 

stagnation point. Since U=0 at x=0, has infinite 
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values for values of z other than those for which 
IV z= X is one of the roots of the numerator of equation 
(11). At x=0, C7=0, the numerator, a cubic in X, has 
roots X! = 7.052, X 2 = 17.80, X 3 = -72.26. The corre- 
sponding values of z are 2i = 0.89, £ 2 =2.25, z z = — 9.13. 

dz 

At z 2 , and z 3 we find singular points at which ^ 

1- indeterminate. A solution curve can then leave the 

2- axis only at one of the singular points. The singular 
point at z 3) involving an imaginary 8, is obviously of 
no interest. It was found that a solution curve leaving 
the singular point at z 2 led immediately to very large 
values of 5, entirely out of agreement with experiment. 



The zero isocline then leaves with zero slope. A little 
consideration will show that the solution curve must 
leave the singular point along an isocline whose initial 
slope equals the value marked on the isocline; in this 
case the one with zero slope. For the sake of accuracy, 
the solution curve was determined from a large dia- 
gram containing more values of the slope than that 
shown in figure 14. The values of z and x resulting 
from the solution are given in table II. From the 
values of V and R, X and 5 were calculated. These 
quantities appear in the same table. The speed 
distributions were calculated from these values of X 
and z. 




ion 



Figure 12.— Speed distribution in the separated boundary layer and wake (experimental). The * indicates an uncertainty in the direction of the speed in the wake. 

Vb means distances normal to the major axis of the ellipse. 



The singular point at z x was the only one that gave a 
solution curve in agreement with experiment. The 
isoclines leaving Z\ have an initial slope given by 2 



2 0.8[ -9072+ 1670.4X- (^47.4+4.8^^) X'- (l+^r) * 3 ] 
L r (-213.12+5.76X+X«) 

a constant for any one isocline. Differentiating with respect to x and setting 

'TO 



(T 



=0,- 



dx 



= 0, U= 0 for x = 0; t hen SU bst i tu tingXi = ( I7')o *i, we gel I he resu 1 1 
t = _0 I80 (^) where (jfi) initial slope of the isocline whose value 

"(g),- 



V. DISCUSSION OF RESULTS 

It will be observed in table II that the boundary 
layer has an initial thickness of about 0.6 millimeter. 
It is a matter of conjecture as to whether a true bound- 
ary layer really exists at the stagnation point. The 
measurements indicate that potential flow probably 
extends to the surface. 

Also of note in table II is the absence of any X less 
than —5.37. The value X= — 12 required for sep- 
aration was not attained. No calculations were made 
for values of x greater than 2.937, but it is apparent 
from the value of U' that no z can be reached for which 
\=U'z= — 12. Therefore, according to the calcu- 
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lation, separation does not occur anywhere on the 
cylinder. As we have seen, separation was found 
experimentally at x=1.99. 

The failure of X to reach sufficiently high negative 
values is reflected in the growing discrepancy between 
calculated and observed distributions in figures 8 to 
11. As shown by figures 4 to 7, Pohlhausen's method 
yields good agreement with experiment up to an x of 
1.832. 

From an estimate of the pressure distribution re- 
quired to give a calculated value of separation in 
agreement with experiment, it was at first suspected 
that the pressure distribution might not have been 
determined accurately enough in the region where 
separation occurred. A redetermination with pres- 
sure orifices installed in the cylinder at critical points 
showed that the pressure distribution as originally 
determined was correct and that as far as 2=2.1 the 

— x b Scale > 




Separated 
^boundary 
' layer 



Reference 
length (L) = 
3.98 in. 
100.9 mm 



6 at staqnotion point from isocline calculation = 0.00602 
Figure 13.— Diagram showing elliptic cylinder and boundary layer in section. 

This diagram shows the thickness of the boundary layer compared to the size 

of the ellipse. 

experimental curve of figure 2 is definitely fixed both 
as regards position and shape. 

The question arises as to whether Pohlhairsen's 
solution begins to fail after x= 1.832 or whether the 
assumptions made in simplifying the general equation 
of section II cease to be valid. The following table 
shows that 8 changes but little through the interval 
in x where discrepancies begin to appear. The failure 
can then hardly be attributed to the growth of <5. 
J. J. Green (reference 8) has shown for the case of a 
circular cylinder that the neglected terms become 
important only at the separation point. Here v 
probably reaches an order of magnitude near that of u. 



X 


b 

(by experi- 
ment) 


0. 1800 


0. 009 


.357 


.014 


.545 


.019 


.725 


.023 


1.097 


.032 


1.457 


.040 


1.832 


.044 


1.946 


.045 


2. 029 


.049 



It appears, therefore, that the initial failure is the 
result of the approximations introduced in Pohlhau- 
sen's solution; whereas farther downstream, say at or 
near the separation point, the failure may be due to 
the absence of important terms in the basic equations 
as well. 



The only other known test of Pohlhausen's solution 
is that given by Pohlhausen himself (reference 1) 
using data obtained by Hiemenz for a circular cylinder 
in water (reference 9). In this case the calculated 
separation point was in good agreement with experi- 
ment. An examination of the results of Hiemenz 
shows a steep pressure rise beginning a short distance 
ahead and extending up to separation. The x, X 
curve for the circular cylinder shows a steep fall in the 
region of separation. Therefore the experimental sep- 
aration point must be determined very accurately to 
show that it falls near \= — 12. Since the foregoing 
condition does not obtain to so great a degree in the 
present experiment, the test applied here to the solu- 
tion in the region of separation is a more sensitive one 
than that applied by Pohlhausen. 
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Figure 14 — Isocline diagram and solution curve of equation 11. 

CONCLUSION 

The experimental data presented here show that 
Pohlhausen's solution may be expected to yield reli- 
able results where the speed in the potential flow is 
increasing. In the region of decreasing speed, the 
solution cannot be relied upon to determine the speed 
distribution, but gives the boundary-layer thickness 
to a fair degree of approximation. Separation may 
actually occur where the solution of Pohlhausen fails 
to give it. 

The author wishes to acknowledge his indebtedness 
to Dr. H. L. Dryden of the National Bureau of Stand- 
ards for his many valuable suggestions and to Profes- 
sors K. F. Herzfeld and F. D. Murnaghan of The 
Johns Hopkins University for their critical reading of 
the manuscript. The author also wishes to express his 
appreciation of the assistance given by Mr, W. C. 
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Table I.— DATA DERI V ED PROM PRESSURE DISTRI- 
BUTION FOR CARRYING OUT SOLUTION OF POHL- 
HAUSEN 



Orifice 
number 



10 



0 

.175 
.260 
.357 
.460 
.545 
.725 
.950 
1.097 
1.250 
1.350 
1.457 
1.581 
1.600 
1.700 
1.704 

1 . s:v2 
1.900 
1.957 
2.000 

2. 079 
2. 100 
2. 196 
2. 56S 

2. 937 

3. 307 



P-Po 



1.000 
. 163 
-. 115 
-.307 
-.440 
-.513 
-.605 
-.657 
-.672 
-.677 
-.677 
-.671 
-.655 
-.650 
-. 627 
-. 626 
-. 592 
-.569 
-.553 
-.540 
-.523 
-. 516 
-.501 
-.497 
-. 509 
-.495 



dp 
dx 



0 

-4. 43 
-2. 56 
-1.571 
-1.010 
-.719 
-.354 
-. 150 
-.0633 
-.0050 
.030 



. 180 
.24 



.306 
.30 



.20 
. 15 
.064 
.001 

.08 



dx* 



-125.4 
26.7 
15.9 
7. 80 
4.30 
2. 72 
1.38 
.60 
.485 
.44 
.47 
.49 



.52 
.50 



0 

-.30 



-.60 



.70 

.70 

.20 

.2 

.2 



U 



0 

.914 
1.055 
1. 142 
1.200 
L 230 
1.266 
1.286 
1.293 
1.295 
1. 295 
1.292 



1. 2S1 
1.275 



L 261 
1.252 



1.240 



I. 230 
1.225 
1. 223 
1. 22N 
1. 222 



7. 92 
2.42 
L 213 
.688 
.421 
.292 
.140 
.0583 
.0245 
. 00193 
-.0116 
-. 0332 



•. 0701 
-.0941 



.121 
. 120 



. 101 



. 0813 
. 0592 
.026 
.000-11 
.03 



UU" 
U'i 



0 

-3.28 
-6. 42 
-9.27 
-13. 14 
-16.9 
-36.4 
-89.4 
-406.0 
-59000. 0 
-1750.0 
-226.0 



-54.0 
-29. 26 



-1.0 
9. 45 



28.6 



52. 1 
97.0 
145.0 
-7. 1X10* 
-134.0 



(a) (7=Vl-(P— Po). 
1 dp 
Wdx' 



(6) U'= 



(c) 



IL" 



Table II.— VALUES OF X AND 5 CALCULATED BY 
POHLHAUSEN'S SOLUTION AT POSITION x AND 
FOR REYNOLDS NUMBER R 

IBoundary-layer thickness in millimeters = 101.1 5| 



X 


V 


z 


X 


R 


5 


0 


7. 92 


0. 890 


7. 052 


24, 500 


0. 00602 


. 180 


2. 42 


1.91 


4. 63 


24, 400 


.00884 


.357 


.688 


4. 22 


2. 90 


23, 500 


.0134 


.545 


.292 


7.26 


2. 12 


24, 000 


.0174 


.725 


.140 


10.47 


1.46 


23, 600 


.0211 


1.097 


.0245 


18.36 


.449 


22, 700 


.0284 


1.457 


-.0332 


27. 48 


-.912 


22, 700 


.0348 


1.832 


-. 121 


42.0 


-5. 10 


24, 300 


• .0415 


1.946 


-. HI 


48.5 


-5. 37 


23. 900 


.0447 


2. 029 


-.095 


52.9 


-5. 03 


23, 600 


.0474 


2. 133 


-.075 


58.5 


-4. 39 


23,400 


.0500 


2. 196 


-. 059 


61.7 


-3. 65 


23, 500 


.0512 


2. 568 


.26 


80.3 


2. 10 


2?., 500 


.0584 


2. 937 


.00041 


• 93.9 


.038 


23, 100 


.0638 



X= U'z 5 = 



O 




Positive directions of axes and angles (forces and moments) are shown by arrows 



Axis 


Force 
(parallel 
to axis) 
symbol 


Moment about axis 


Angle 


Velocities 


Designation 


Sym- 
bol 


- !^ 
Designation 


T, ' -/. .;. 
Sym- 
bol 


Positive 
direction 


Designa- 
tion 


Sym- 
bol 


Linear 
(compo- 
nent along 
axis) 


Angular 


Longitudinal 

Lateral 

Normal 


X 
Y 
Z 


X 
Y 
Z 


Rolling 


L 
M 
N 


Y >Z 

Z >X 

x — >r 


Roll 


e 

* 


u 

V 

w 


V 
Q 
r 


Pitching 

Yawing 


Pitch 

Yaw 











Absolute coefficients of moment 
(rolling) (pitching) 



Cn ~qbS 
(yawing) 



Angle of set of control surface (relative to neutral 
position), 5. (Indicate surface by proper subscript.) 



4. PROPELLER SYMBOLS 



D, Diameter 

p, Geometric pitch 

p/D, Pitch ratio 

V', Inflow velocity 

V 8 , Slipstream velocity 

T, Thrust, absolute coefficient C T = 

Q, Torque, absolute coefficient C Q = 



1 hp. = 76.04 kg-m/s = 550 ft-lb./sec. 
1 metric horsepower = 1.0132 hp. 
1 m.p.h. = 0.4470 m.p.s. 
1 m.p.s. = 2.2369 m.p.h. 



a, 

V, 

n, 



T 

P n 2 Lh 

Q 

pn 2 D 6 

5. NUMERICAL RELATIONS 



Power, absolute coefficient = pn zjys 

Speed-power coefficient = -^p^ 
Efficiency 

Revolutions per second, r.p.s. 
Effective helix angle = tan*" 1 Q^;) 



lib. = 0.4536 kg. 

1 kg = 2.2046 lb. 

1 mi. = 1,609.35 m = 5,280 ft. 

1 m = 3.2808 ft. 



